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On the speed of convergence in the local limit 
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Abstract 



^^ ■ We establish the upper bound on the speed of convergence to the infinitely 

divisible limit density in the local limit theorem for triangular arrays of random 
Q_il variables {Xk^n, k = 1, ..,an, n € IN}. 
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Ph 1 Introduction 

C^ I This paper is motivated by [KKu07], where the local limit theorem for a triangular 

array of independent and identically distributed (i.i.d) in each series random variables 
{Xkni k = 1, .., ttn} is established. Staying in frames of the situation studied in [KKu07], 
in this note we would like make a step further and obtain the information about the speed 
of convergence to the limit density. 



0\ ' In contrast to the local limit theorem for the normal law, there is not much known even 



about the local limit theorem for infinitely divisible limit densities. Of course, one can 
refer to Gnedenko's theorem on the necessary and sufficient conditions for the convergence 
^ ! to the stable law, see [IL71]. Under certain conditions the uniform convergence to the 

^^ I limit density was proved in [KKu07], but up to the author's knowledge in the general 

case nothing is known about the speed of convergence. 

To make the presentation self-contained, we quote below the necessary and sufficient 
rS • conditions for convergence to the infinitely divisible law, see Theorem 2, Chapter XVII 

C^ : §2 from [Fe71]. 

Recall that a measure M on R is called canonical if M(I) < oo for any finite interval, 
and 

M+(x) = / —Midu) < +00, M-ix) = / —Midy) < +oo, x > 0. 

Jx W J-oo U^ 

A sequence of canonical measures {M„} converges to a canonical measure properly if 
Mn{I) -^ M{I) for any finite interval, and M+(x) -)■ M+(x), M~{x) -)■ M~{x) for every 
X > 0. In this case we write M„ — )■ M. 
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Theorem 1. [Fell] Let {Xk.n 1 < ^ < o-n] be such that Xk^n o-f^ i.i.d. for any 1 < k < an, 

a„ — )■ oo as ra — 7- oo, and satisfy 

limP{|Xi,„|>£} = 0. (1.1) 

n— >oo 

Let Fn{du) he the distribution function of Xi^n, 

Mn{du) := anU^Fn{du), l3n '■= / sin uFn{du). 

Jn 

Then 5'„ := Xi „ + ... + Xa^^n converges in distribution to some random variable S if and 
only of 

Mn -> M, a„/3„ -)■ /3, as n -)■ oo (1.2) 

for some /3 G R and some canonical measure M . In this case the characteristic function 
^{z) of S is given by 

$(^) = exp {i(3z+ / M{du) } . (1.3) 



R 



U 



The function 



^{z) := -iPz + / X M{du) = -tPz + (P{z) (1.4) 



R U^ 



is caUed the characteristic exponent of the infinitely divisible variable S. Put 

Ipniz) := -2/3n2: + / 5 Mn{du) = -if3nZ + 0„,(2). (1.5) 



Remark 1. Of course, one can formulate Theorem 1 with the function 1|m|<i instead of 
sinw under the integral, but for technical convenience we need the Levy representation 
(1.3). 

Sometimes, especially when the convergence in Theorem 1 is that to a stable law (cf. 
[IL71]) it is more convenient to consider the random variables in the form 

X,,„ = 1^, (1.6) 

where the variables ^k,n, ^ ^ k < an, are i.i.d. for each n, and the sequence (6n)n>i 
satisfies certain growth assumptions. In the paper we assume that the random variables 
Xk^n are of the form (1.6). We assume that the conditions of Theorem 1 hold true, and 
thus Sn = f^' converges weakly as n — )■ oo to some infinitely divisible distribution 

S. Under some conditions on the sequences (a„)„>i, (fc„)n>i, and on the distribution of 
.^i^„ (cf. [KKu07]), Sn and 5* possess transition probability densities and the local limit 
theorem takes place. Taking this result as the starting point we derive in Theorem 2 
(under some additional assumptions) the speed of convergence to the limit density, and 
illustrated our result by examples. 

In order to make the presentation as transparent as possible, we write the main 
notation in Table 1. Finally, denote by m the symmetrization of the measure m, i.e. 

^(^) ._ rn(A)+rn{~A) ^^^ ^^^ g^^^^ ^^^ A G R. 



Table 1: Notation 



Variable 


Char, function 


Probab. measure 


Probab. density 


Sn 
S 


0niz) 
Oniz/bn) 

$(z) = e-'^(^) 


Gnidu) 

Fn{du) = Gnipndu) 
Pnidx) 

P{dx) 


gn{.u) 

fn{u) = hngn{hnU) 
Vn{x) 

p{x) 



2 Main result 



We assume that the assumptions below hold true: 

A. for any n > 1 the variable ^i „ possesses the density gn{x)', 

B. 3a G (0,2) such that Keipiz) > c\z\"' for \z\ large enough; 
C.y6 > Owe have N{6):= sup |^„(2)| < 1; 

">l,|z|>'5 

D. sup„>i J^gl{x)dx < oo; 

E. bn -)■ oo, ^sAl _i. as n -)■ oo; 

F. for n > 1 one of the conditions below is satisfied: 
a) there exists c{6) > 0, < k < 2, such that 

On / (1 — cos{zu))Fn{du) > c{6)\z\'^ for all \z\ < Shn] 
Jn 



(2.1) 



h) F < Fn on R. 

G.35 >0 such that inf iReOnizM > 0. 

n>l,\z\<S 

Remark 2. Condition A, and C-E are taken from [KKu07]. Instead of condition B in 
[KKu07] another condition was assumed (namely, a version of the Kallenberg condition 
[Ka81] for the sequence of measures M„), which in fact implies B. 



Let 



In ■ = sup 



\Re(f){z) — Re(f)n{z) 



In 



|Im0(2;) — Imcpniz) 



f2.2l 



l + z^ ' '"■ .m l + z^ 

\n • \(^nPn P|; 

where /3„ has the same meaning as in Theorem 1. 

From now we fix 5 > 0, for which the above conditions hold true. For some fixed 
< e < 1 put 

p,4n) := max (x„, 7^, In, a-\ e'^"(i-^W+^)e-(i-)^^'^(^^")) , (2.3) 

where N{6) is defined in C. 

Theorem 2. Suppose that conditions (1.1), (1-2), and AG are satisfied. Then the 
distributions Sn and S possess, respectively, the densities Pn{x) andp{x), and 



sup \pn{x) -p{x)\ <Gp^^s{n), n^oo, 



(2.4) 



where pe,&{n) is given by (2.3). 



One can simplify the expression for the speed of convergence, but at the expense 
of some additional assumptions on a„ and 6„,. We say that a sequence (c„)„,>i satisfies 
condition H, if there exist a,b > such that 

< liminf — < limsup -^ < oo. 

Corollary 1. Suppose conditions of Theorem 2 hold true, and assume in addition that 
the sequences {an)n>i o-nd {bn)n>i satisfy H. Then 

sup|p„(x) -p{x)\ < Cmax (7n:7n, Xn, <M- (2.5) 

Corollary 2. Suppose that conditions A F and H hold true, the densities Pn{x) and 
p{x) are symmetric, and 

^niz) > $(z) Vn > 1, (2.6) 

uniformly in {z : \z\ < 6bn}- Then 

sup \pn{x) - p{x)\ < C(7^ + r(n))), (2.7) 

where r{n) = o{n~^) as n ^ oo for any k > 1. 

Remark 3. As one can expect, the oscillation of measures involved in 7^ and 7^' can 
play the crucial role in the estimation of the speed of convergence. For example, it might 
be insufficient to know the behaviour of such a "rough estimate" for 7^ as below: 

/j^(l A \uz\'^){Mn{du) - M{du)) 



sup 

2GR 1 + Z^ 

in particular, when the densities {gn)n>i have oscillations. Such a situation is illustrated 
in Example 1. 

3 Proofs 

Proof of Theorem 2. Recall that the densities p.„(x) and p{x) can be written as the inverse 
Fourier transforms of the respective characteristic functions: 

p{x) = (27r)^i f e~'''=^{z)dz = (27r)-i f e-''''-^^'^dz, (3.1) 

Ju Ju 

p„(x) = (27r)-i / e-''^<l>rriz)dz. (3.2) 

By (3.1) and (3.2) we have 

An: = 2TTsnp\pnix)-p{x)\< / l^niz) - ^{z)\dz 
xeR Jr 

< ( /" " l^niz) - Hz)\dz + f \^niz)\dz + f |<l>(z)|dz 

J ~5b„ J\z\>Sb„ J\z\>Sb„ 



~Sbn J\z\>5bn J \z\>5bn 

--: Un) + hin) + h(n) 



where S > 0. We estimate the terms Ik{n), k = 1,2, 3, separately. 
Estimation of Ii. Observe, that 

|1 - e^'+'^l = |1 + e^^' - 2e^' cos 1/1^/2 

= 1(1 -e")2 + 2e"(l- cosy) 1^/2 

< |(l-e")2 + eVr/' 

< |1 - e""! +e^^'^\y\ 
<e^^{\x\ + \y\), 

where x, ?/ G R, and x+ := max(x,0). Denote 

z 
Hn{z) := ?/;(2;) + a„ln6'„(— ). 



Then by (3.3) we get 



rSb„ 



/a On 
g-ReV.(.)+(ReH„(.))+(|J^g^^(^)| + |Inii7„(z)|)ci2. 
'Sbn 

Since ln(l — z) < —2; for z G (0, 1), then 

ReHn{z) = Retlj{z) + a„ln |^n(f )| < Re^(z) - a„(l - |^n(f )|). 



Observe, that 



6'n(7-) = / cos{zu)Fn{du) = a^^ / 

"n JR Jr 



COSUM 



Mn{du). 



'R JR ■" 

Therefore, by (3.5) and F we have for all n large enough and \z\ < 5hn 

Re^(^) _ [ReH^{z))^ > a„(l - |^n(f )|) 



= a„ / (1 — cos{uz))Fn{du) 
Jn 

>ci5)\zr. 



if F.a) holds true, or 



Rei/jiz) - {ReHn{z))^ = Reip{z), 
if F.b) is satisfied. On the other hand, for z G (0, 1) we have 



iu^+i-.i<5:<i^<i^Bi--)' 

fc=2 k=l 



< 



2z 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



Then by (3.6) and G we derive 
\ReHn{z)\ < ReV'(^)-a„(l-|^„(^/6„)|) 



{l-\en{z/br,)\)-\n\en{z/bn)\ 



< 



[1 - cos{zu))u~^{M{du) - Mn{du)) 



R 



+ 2 ^an( / [1 — cos{zu))Fn{du)\ ■ { I cos{zu)Fn{du) . 

^Jr ^ ^i]R ^ (3.9) 

< ci (7^(1 + z^) + {2an)-^{ReM^)T) 

< ci (-fUl + z^) + (2a„r 1 ((1 + z')^n + (1 + ^') 

Next we estimate |Iniif„(z)|. Observe that for z = x + iy, where x, y G R, 

Inline = Arg z = arctan— , 

X 

and for all x G R we have | arctan x — x\ < C3|x|^, where C3 > is some constant. 
Therefore, 

\lmHn{z)\ = Imipiz) + a„Im ln6n{z/bn) 



< 



< 



+ 



i3z 



2;sinM — sin(2M) lm.9n(z/bn 

2sinM — sinf^w' 



Pz + z -^M^{du) 
'r ■" 



+ 



R 



R 



smi zu) — zsmu 



sin u , ^ , , , 
z^-Mn{du) 



Mn{du) 

1 



Re9n{z/bn 



Re6niz/br^ 



--1 

arctan 



{Mn-M){du) 



lm9n{z/bn) lm9n{z/bn) 



Re9niz/bn) Re6n{z/br, 



< hi + /12 + /l3 + hi + /l5- 



Observe that by M„ — )■ M and a„/3„ — )■ /3 (cf. (1.2)) we have 

hi{n) < \z\\/3 - an/3n\ = \z\Xn- 

For /i2(^) we have 

/i2(n) < |lm0,(z) -Im0(z) I < 04(1 + ^^)7;' 

Using G, we derive 

hsin) < cga;;;^ |Im 0„ (2) I Re 0„ (2;) < C6a~'^{l + z^f. 

Analogously, 

/i4(n) < CTl3na~'^\z\Re(j)n{z) < csa;^'^\z\{l + z^). 



(3.10) 

(3.11) 

(3.12) 

(3.13) 
(3.14) 
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Finally, for J15 we have 

/l5 < a„C3 



< CgOn' 



Cga^' 



lm9n{z/bn) 



< cga^ 



Re9niz/bn 

sm{zu) — 2 sin It 



sin{zu) Fn{du) 



R 



Mn{du) + z 



smu 



R 



n^ 



Mn{du) 



luKpniz) - Zf3„ 



<c,oa-\l + z'f. 
Thus, we arrive at 



hin) < ciimax(fi;„,7^,7^',a„ ^ 

< C12 max(K„, 7^, 7^', a'^). 
Estimation of I2. We have by C and D 



Sbn 



<^>\l + z'^fdz 



(3.15) 



h(n) 



2|><56„ 



Or. 



dz = br, I \6n{x)\''"dx 

\x\>5 



<bnN{6r"-' / \eM\'dx 

J\x\>5 



< bnN{5Y--' sup / gi{x)dx 

n>l Jia 

Take e > such that In A^(5) + £ < 0. By E, ^^ -;■ as n -> 00, and thus without loss 
of generality we may assume that for all ri > 1 we have ^^Al < ^_ Then 

^ ga„ln7V(<5) ^ g-a„| In Ar(5)+e| 

Estimation of I^. For any e > 

-(l-e)Rei/>{5fe„) 



h{n) < cu [ e-'^'^^^'Uz < ci5(e)e- 
Summarizing the estimates for Ii{n), i = 1,2,3, we derive A,„ < Cp^^s{n). 



U 



Proof of Corollaries 1 and 2. Clearly, the proofs are obtained as slight modifications of 
the proof of Theorem 2. Since a„ and 6„ satisfy H, the terms /2(^) and I^in) decay as 
o(ra~^), n — 7- 00, for any /c > 0. This implies the statement of Corollary 1. To complete 
the proof of Corollary 2, we need to estimate more precisely Ii{n). Let us look closely on 
the properties of the function Hn{z) from (3.4). Since both p„(x) and p{x) are symmetric, 
Hn{z) is real-valued. Further, condition (2.6) implies that Hn{z) > 0. Therefore, instead 
of (3.9) we get 

Hn{z) < ij{z) + anlnOniz/bn) < 7n(l + Z^), 

which implies 

h{n) < a^. 

n 



4 Examples 

Example 1. Let (^n)n>i be i.i.d. random variables with probability density 

(1 — cosm) 
q(u) = Ca , M, -, u eR, < q; < 2. 

Then one can check (using Theorem 1 with an = n and 6„ = n^/°) that 

Q Si "T •• "T Sn Q 

'^'^ •- :^ ^ '^' 

where S* is a symmetric a-stable distribution. In this case the respective measure M{du) 
in (1.3) is equal to Ca\u\^~"du, and after the appropriate choice of Ca we have ip{z) = |2|". 
For example, in the case a = 1 we must chose Cq = l/ir. Clearly, conditions A, B, D and 
G are satisfied. Condition C is the Cramer condition (cf. [Lu79]) for the characteristic 
function of ^i, which is satisfied since the law of ^i is absolutely continuous. 
Let us check condition F. Consider 

"^l-cosizu) ^, f°° 1-cosu 1- 

cos{n ' u)du = \z\ / — cos(ra ' u/z)au. 



" Jo ^ 



We need to estimate from above 



I{a,k):= I j— — cos{kv)dv. 







V 



Note that for a = 1 we have (cf. [BE69], p. 28) 

/(l,^) = f(l-|^l)+- (4.1) 

It is also possible to calculate I{a, k) for a G (0, 2)\{1}. Integrating by parts, we get for 
any A; > 



sin kv 1 — cos v 



°° 1 f°° ■ . wsinw - (1 + a)(l - cost;) 



/(a, k) = —r -— — ^ - T / sm{kv) -— dv 



k In V 



k V 

1 f°° ., , tisinf — (1 + a)(l — cosf ) , 

■- / sm(kv) — dv. 

kJo v^+^ 



^ , , , , smikv) sinf , 
h{a,k):= -— dv 



2+a 



The integrals 



and 

hia,k):=J^ ^^ dv. 

can be calculated explicitly, see [BE69], p. 77- 78, from where one can derive the asymp- 
totic behaviour as /c — )■ oo: 

h{a, k) = j l,,:,.\ ,' - c^Akr\ (4.2) 

4 r(l + a) sm(7ra/2) 
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hia, k) = 2"2r(-l - a) cos(7ra/2) [2|A;|°+^ -\k+ 1|°+^ -\k- 1|"+^] ~ Ca,2\k\''~\ (4.3) 

where c„,i := 2r(i+a) sm(7r«/2) ' ^".2 = 2~^a(« + l)r(-l - a) cos(7ra/2). Thus, we have the 
exact expression for I{a, k), from which we derive 

/(a, k) ~ Cs^Q./c"""'^, /c — !■ 00. 

where cs^q, = ci^q, + 2(a + 1)c2,q. Finally, for \z\ < 5n^^°' with 5 > is small enough, 

r - ( ■n}^'^ \ 

an / {I - cos{zu))Fn{du) = |^|° l-2cj{a,——) > c{6)\z\'', a G (0,2), 

where c{6) > is some constant. 

Let us calculate the order of convergence. From above, we have for a G (0, 2) 



Jn m 
Thus, by Corollary 1 we arrive at 



Cz 
< 



^{2-a)/a ■ 



Pin) < L, _^ . ^ ^ ; (4.4) 

\Cn a , 1 < a < 2. 

Example 2. Suppose now that ,^i^„ possesses the distribution density 

9n{U) := 2nusinh(«/n)-'^l"l>l' 

and ttn = bn = n. Conditions A, C-E were already checked in [KKu07], in particular, 
it was shown that Sn converges in distribution to a hyperbolic cosine distribution S, i.e. 
the distribution density of S is p(x) = — ^— . Since in this case 

•' ^ ^ ' n cosn X 

anfn{u) = - — r-T— l|„|>i t 7^ — r-7— =: f{u) as n ^ oo, 
2nsmhn ' '-« 2nsmhn 

the function 

■iIj{z) = / (1 — cos{uz))f{u)du 
Jn 

satisfies condition B with a = 1. Let us check F for k = 1. Since for \z\ < 1 we have 
1 — cos z > {1 — cos 1)2;^, then estimating ^.^^ from below for small n by a constant we 
get 

a„ / (1 — cos(zn))/„(u)(iu > 72(1 — cos 1) / {zuf fn{u)du 

Jr J\uz\<1 

> fl — COS llkl inf \z\ I — : — du 



\<5n ./i/n sinhu 



> Ci\z\ mr U 1— r 

\z\<Sn \\z\ n 

>ci(l- (5)1^1, 



uniformly in {z : \z\ < 6n}. Thus, condition F holds true. 

It remains to check condition G. Let \z\ < d. Since the function rsinh(u/r) is in- 
creasing in r, we have by dominated convergence theorem 



inf l^nf^)! = inf 

n>l n>l 



COS{ZUj 

„>i riM sinh(M/n) 



°° cos(zu) 



lim \Oniz)\ = I — du 



> 1 r^'^" 

which gives G. 

Finally, by Corollary 1 we arrive at 

C 
sup|p„(x) -p{x)\ < —. 

Acknowledgement. The author thanks A.Kulik, who initiated this work, for helpful 
discussions and comments. 
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